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Abstract 

We investigate the existence of bound states of baryons in a kaon condensate 
using chiral mean field theory. The interactions are described by an effective 
SU(3)lxSU(3)p{^ chiral lagrangian where terms of higher order in density, baryon 
momentum, and kaon mass are suppressed by powers of the symmetry break- 
ing scale, A. We take up to next to leading order terms (n = 2,3,4). We 
search for infinite baryon number solutions, namely "strange baryonic matter", 
using a Thomas-Fermi approximation for a slowly varying condensate and a 
lowest order Hartree approximation to describe the many body interactions. 
For simplicity we study a pure condensate and only neutrons, the lightest 
baryons in that condensate. We find solutions with neutron number densities, 
Pn ^ 3.5/30) where po is the infinite nuclear matter density. This is consistent 
with the estimate of the onset of a K-condensate at pn — 2-4 pQ. We show that 
the binding energies, Eb, grow with pn and for p„ < 7po (at pn ^ 7po pertur- 
bative expansion is lost) we find < IbOMeV {E^ < JOMeV for pn < 5po) 
even in the most favorable cases. These binding energies may be too low for 
this type of matter to appear and persist in the early universe. 
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1. Introduction 



QCD is the theory of strong interactions among quarks and gluons, the elementary 
particles that constitute hadrons. With three massless quarks, m„ = = = 0, 
the QCD lagrangian has a global chiral SU(3)LxSU(3)j^symmetry. This is considered 
to be an accidental symmetry, since there is no deeper reason for the (almost) mass- 
lessness of these three quarks. The chiral symmetry is approximate since the three 
quarks masses are small, but not zero. The symmetry is most explicitly violated 
by the s-quark mass, ~ 24:0MeV, which is much larger than the u and d masses 
(rriu ~ 6 MeV, rud — 12 MeV). Even if we can not derive from QCD the properties of 
the quark and gluon bound states, the hadrons, we know the mass spectrum and low 
momentum interactions of hadrons show the persistence of the SU(3)lX SU(3)j^ chiral 
symmetry. The SU{3)v vectorial and SU{3)a axial vector symmetry groups (those 
whose generators are the sum and the difference respectively of the left and right 
generators) are realized differently. The vectorial subgroup SU{3)v is realized in the 
Wigner-Weyl mode, yielding (almost) degenerate multiplets (octets and tenplets) of 
baryons and mesons. This is the (approximate) symmetry used to classify hadrons 
(in the "eightfold way"), that lead to the proposal of quarks as a means to populate 
the fundamental representation of the group. The axial symmetry SU{3)a is instead 
realized in the Nambu-Goldstone mode, namely it is spontaneously broken at a scale 
A ~ IGeV yielding an octet of (quasi) Goldstone bosons, one for each broken gener- 
ator, the pions, kaons, and eta mesons. The lightness of these mesons compared to 
the other hadrons justifies this identification. The mass of the Goldstone bosons is a 
result of the explicit breaking of SU (3) a due to the non-zero quark masses. In fact 
from the lowest order (n = 2, see below) chiral lagrangian, the tt, K, and rj masses 
result linearly proportional to the u,d, and s quark masses. At this order one gets the 
phenomenologically successful Gell-Mann Okubo relation among meson masses. This 
is evidence that the perturbative expansion in the chiral lagrangian is good, namely 
that higher order terms (that modify the Gell-Mann Okubo relation) are small. 

A non-linear effective chiral lagrangian |I| is the most general lagrangian for 
the baryons and the octet of quasi Goldstone bosons (therefore valid at energies 
below the scale of spontaneous chiral symmetry breaking, i.e. E < A) which is 
compatible with the approximate accidental SU(3)LxSU(3)j^chiral symmetry of QCD 

H. Although there are different ways of parametrizing the Goldstone bosons, it 
has been shown that they all lead to the same observables This effective non- 
renormalizable lagrangian consists of a power series expansion in derivatives, baryon 
fields and the chiral symmetry quark mass matrix A4 = diag{my_, rud, rUs)- As a trick 
to use chiral symmetry to also fix the form of the explicit symmetry breaking 
terms using our knowledge of QCD, the matrix Ai is promoted to a field with its 
chiral transformation chosen so that it would fix the form of the Ai dependent terms 
in the QCD lagrangian (if chiral symmetry were exact). The chiral symmetry breaking 
terms in the effective lagrangian lead to s-wave interactions of the Goldstone mesons. 
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that are very important in the phenomenon of kaon condensation ^. Without 
exphcit breaking the Goldstone bosons have only derivative couphngs, which can be 
understood by recalhng that the Goldstone fields are the angular coordinates that 
parametrize the orbits of degenerate vacua, so that a change in value of the field can 
not affect the energy. 

A chiral effective lagrangian includes all the terms compatible with the approxi- 
mate chiral symmetry of QCD with coefficients to be determined phenomenologically 
when possible. Terms with dimension larger than four have dimensionful parameters 
which are proportional to, either, inverse powers of the symmetry breaking scale, 
A ~ IGeV, or inverse powers of as shown in Eq.( |2.1| ) below. This prescription in- 
sures that loop corrections to each term generate terms of the same form, if A ~ Airf.,^ 
0, Q . The usefulness of the expansion resides in the ability of cutting the series after 
a few terms. Assuming that the expansion parameters (which are, d^/A, rris/A, and 
Pb/ f^A = pb/7po, see below) are small, the terms in the lagrangian can be organized 
in successively less important sets and a finite number of terms can be used. 

This method provides the only systematic way to implement the symmetries of 
QCD in TT-i^'-T^-baryon interactions. This method has proven to be useful, in many 
applications, among which describing properties of bulk hadronic matter, such as the 
formation of a kaon condensate [§, 0. It has been shown that due to the kaon's 
dominant s-wave coupling to baryons the formation of a kaon condensate is quite 
insensitive to nuclear interactions. It is believed that a kaon condensate will most 
likely form at a baryonic density anywhere from 2-4 times nuclear density. 

In this paper we are not interested in the details of the onset of the condensation 
but rather in solutions to the non-linear classical field equations that describe an 
isolated system consisting of baryons in a K-condensate. In order to approach such a 
complicated problem many approximations have to be made. The result is a system 
of equations that resembles those of a liquid droplet ^ where the baryons are 
treated as a gas of pseudo particles trapped in the bose condensate of kaons [|10 



Lynn, Nelson, and Tetradis ||T0| (LNT from now on) studied this problem using a 
phenomenological combination of a chiral lagrangian, having no terms with four or 
more baryon fields and no terms which are higher order than linear in quark masses, 
and the Walecka lagrangian. The Walecka model |]ll| consists of two fictitious massive 
vector and scalar fields, and 0, coupled to protons and neutrons in a renormalizable 
lagrangian. It describes well the properties of bulk nuclear matter. LNT added 
the Walecka lagrangian to incorporate nuclear forces not included otherwise in their 
model. They also coupled the field to the mesons. Here we only use the chiral 
effective lagrangian with the four baryon terms and next to leading order terms as 
well, that include more than four baryons and higher powers of the quark masses. 
In fact we have shown elsewhere fl^ that the four-fermion terms incorporate into 
the chiral lagrangian the same description of bulk nuclear matter contained in the 
Walecka model, when baryons are restricted to the nucleon doublet. Because the 
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solutions we find have large densities, we need not only investigate baryon-meson 
and baryon-baryon interactions but also three baryon interactions (ie. terms with six 
fermion fields). We are guided here by the belief in a perturbative series dictated 
by the broken chiral symmetry as explained above where the lowest order terms 
are dominant. We think that this is the essential difference between our work and 
LNT's treatment of the same problem, namely we rely on a perturbative expansion 
as described above and, as we argue towards the end of this paper in section 8, we 
believe they do not. 

This paper is organized as follows. The SU(3)lX SU(3)j^chiral lagrangian we use 
is presented in section 2 and the necessary approximations and our ansatz are given 
in section 3. The following four sections discuss the solutions. In section 4 the 
general method to obtain solutions is explained using just the lowest order lagrangian. 
Section 5 shows how the requirement of a continuous density allows us to constrain 
the solutions to bands in binding energy-density space (shown in Fig. 3), with just the 
lowest n = 2 lagrangian. In section 6 and 7 we add higher order terms one at a time 
to see if those terms help to obtain solutions with a lower density for a given binding 
energy. Section 7 contains our main results. In section 8 we compare our work with 
LNT and section 9 contains our conclusions. 

2. The Lagrangian 



A term in the Lagrangian consistent with naive dimensional analysis p| is given by 



Cf!A2 





a 


"n" 


b 




c 


m" 


fM 








A 




.A. 



(2.i; 



where = 93MeV, 



m 



{mx,fnB) and at first order in rris 



m 



K 

24:0MeV, we take m„ = = 0, and F^ represents 1, 75, •y^, 7^75, and a^''' , either 
multiplied by the SUv{S) generators, T^, or not. The order of the terms is given 
by the index a + c + d = n. Terms with n = 2 are unsupressed. Higher orders are 
suppressed by A^""'. 

The meson matrix 11 is defined as HaTa where IIq is the meson octet of Goldstone 
bosons, and = Xa/2, where Xa are the eight Gell-Mann matrices. 
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The baryon octet Ba also appears in the combination B = BaTa, 
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(2.3) 



We consider only the octet of baryons because they are the hghtest of the baryon 
multiplets. Finally, C in Eq.( p.l|) is a coefficient of 0(1). Notice that powers of the 
meson field, 11, are not suppressed by powers of A, since it appears in the combination 
n//7r. The derivative factor operates on both meson and baryon fields. For the baryon 
fields, however, only the spatial derivatives should be included. We are not using a 
heavy baryon formalism here since the effective mass of the baryons will be small 
compared to their momentum (see the examples provided by the values of m* in 
the Table 1). Since we are using a Hartree approximation, where the baryons are 
treated as free pseudo particles, we can use the equation of motion to replace time 
derivatives by spatial derivatives in the interaction lagrangian Notice that the 
density, p = B'joB, appears in the expansion parameter, p/ f^A ~ p/7po where po is 
nuclear density, p = 1.28 lO^MeV^. Therefore, in order to keep the chiral expansion 
reliable we need solutions with densities p < 7po and we do not expect them at p 
below the onset of K condensation estimated at 2-4po- Thus we work in the region 
in which the perturbative expansion in baryonic density is marginal. With p = 3.4po 
the expansion term BTaB/ f^A is of the same order of magnitude as itik/A. 

We choose the terms in our lagrangian which satisfy the condition a + c + d = 2, 3. 
CcMrai = ^f^Trd.J^d^^^^ + ^f^ATrlMJ^ + h.c] 

+ TrB{i^ - mB)B + iTrB-f''[V^, B] 

+ DTrBYldA^ B} + FTrBYl5[A„ B] 

- aiTrB{iMi + h.c.)B - a2TrBB{^M^ + h.c.) 



- C - - 
a^TrlMT. + h.c.]TrBB + -^TrBTo^BBT'^B 

fn 



C H 

+ -^TrBBBV^BBV'B + --^Tr[BV JBBV'B + h.c] 

+ ■■■ (2.4) 

where the dots indicate terms of higher order and a few of the same order, which 
involve different ways of contracting the baryon octet that within our approximations 
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give no new terms later. The non-linear sigma field, S, is given by 

S = exp (2zn/^) (2.5) 

and the field, ^, is defined as 

^ = VE = exp{m/U). (2.6) 

The explicit symmetry breaking is expressed as an expansion in powers of the small 
quark mass m^, an element of the quark mass matrix A4, 



(2.7) 




Finally, we have the meson vector and axial vector currents 



A = '^i^^d,^~m^)- (2.8) 



Most of the coefficients in Eq. ( |2.4| ) are fixed by low energy NN and vrA^ scattering, 
and by mass splittings in the baryon octet. Their tree level values are 

D = 0.81, F = 0.44, 

ai = 0.28, aa = -0.56, = 1.3 ± 0.2, (2.9) 

rus = 240MeV, tub = 584MeV, A ~ IGeV 

and, Ca, Ga, and Ha are free parameters of 0(1). 



3. Ansatz and Approximations 



We choose a simplifying ansatz []T0[ with only one non-vanishing meson expecta- 
tion value, (K^)- In the presence of this VEV the baryon masses are modified with 
the lightest baryon being the neutron. For simplicity we will study the formation of 
bound states with only neutrons. We hope that this can be instructive in search- 
ing for solutions with more degrees of freedom. Our ansatz for the classical meson 
expectation value is, therefore, 




(If) 





n = ^ I 0^ {K^') I (3.1) 
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where {K^) is chosen to be real, and from the hermiticity of 11 we get {K^) = (K^)- 

We are interested in finding classical solutions of the equations of motion consisting 
of bound states at zero pressure and with large baryon number. In order to minimize 
their surface energy these solution should be spherically symmetric. Furthermore, we 
will consider bound states large enough for the surface effects to be negligible. We 
take the density of baryons and the value of the meson field to be almost constant 



throughout the interior. Using the Thomas-Fermi approximation |]I4[ we assume that 
the kaon fields vary slowly compared to the baryon wavelength, i.e. VK'^/fj^ ^ 
k, where k is the neutron's wavenumber. This allows us to effectively treat the 
neutrons at each point as a Fermi gas in a constant kaon field. Finally, since the 
kaon field is slowly varying we find that all p-wave and higher derivative interactions 
are negligible compared to s-wave interactions. With the ansatz given in Eq. (|3.1|) , 



the vector current, V^, vanishes and the terms with axial vector coupling — D and F 
terms — to the neutrons given in Eq.( p.4|) also vanish. With these simplifications the 
lagrangian reduces to 



chiral 



cos( 



(72 

-\-n{i(j) — m*^)n j^n^nn^n 

Cl Gi 
H — —nnnn + _ „. . nnnnnn 



+ -;r:;-^n'nn'nnn + -TT^-rni'y ■ vnnn 



PA 



(3.2) 



m 



K 



moA, n is 



where 6 = ^/2{K^) / f^^, we have used the first order (in m^) relation 
the neutron field, m* is the neutron mass within the (K^) condensate [rU|, 

m* = TTT-AT — afms(l — cos 6^), (3.3) 

with at = 2a3 + 02 + 01 = 2.32 ± 0.4 and mjv the free neutron mass, mjv = 939MeV. 
Notice that a constant has been added to the lagrangian in Eq.( |3.2| ) so that the 6 
dependence of Cchirai disappears from the lagrangian when 6 = 0. This is the origin of 
the (1 — cos6') factor in the second term of Cchirai in Eq.( |3.2| ). Notice that there is also 
a factor (1 — cos^) in m*, Eq.( |3.3|) , so that m-Ar is the physical neutron mass (outside 
the condensate). The constant mjsf contains, therefore, all contributions of the form 
TrBB in Eq. ( p.4|) obtained by setting ^ = 1 and S = 1, thus ttin = 771^ + 27713(03 + 02) 
(actually, knowing 02 and 03 this relation fixes tub, Eq.( |2.9| )). 

The three vector momentum and spin dependent terms containing currents of the 
form and 777*7577, where i are the spatial degrees of freedom, average to zero in 
spherically symmetric bulk matter. Furthermore, all terms containing time derivatives 
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of the condensate vanish as a resuh of minimizing the thermodynamical potential at 
fixed electric charge: the time dependence turns out to be simple harmonic with the 
frequency equal to the electric charge 0. 

4. Infinite Solutions with Lowest Order n = 2 Terms 



Let us first analyze the n = 2 terms of Eq. (pl2|) {Gi = G2 = Hi = 0). 
The Thomas-Fermi approximation allows one to take a free particle wave function, 
expi{k ■ X — tkt), for the neutron, where k and are the space dependent wavenum- 
ber and energy respectively. The Dirac equation for a neutron moving in the mean 
field of all other neutrons is 



a 



CI 



J TV J TT 



n = 0, 



(4.1) 



and by substituting in the neutron wave function, we get 



7°(e, - ^{n^n)) - - {ml - ^(nn)) 



n 



0. 



(4.2) 



Squaring the above equation and applying Dirac algebra we get the dispersion relation 
of a free quasi-particle, 

(4.3) 



,*2 



e--F-m*2 = 0, 
where we define the effective energy and mass of the quasi-particle to be 

* _ t 

fn 



m 



CI 



nn. 



(4.4) 
(4.5) 



The zero temperature ground state, |\1'), of our system can be found by minimizing 
the thermodynamic potential, VL. For fixed baryon number, B, 

Q = Etotai - f^BB, (4.6) 

where /ib is the baryon chemical potential and B = J dP'xn^n. The total energy of 
the system, Etotah is given by 

Etotai = TOO = j ^3^^^^ ^^^^ (4 7) 
where Sk is the kinetic energy density of the condensate, 

= \h'{e' + {vef) = Im^er, (4.8) 
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En is the neutron energy density. 



En 



and Ue is the potential energy density of the condensate 

Ue = f^mlil-cosO). 



(4.9) 



(4.10) 



The thermodynamic potential, Q, is a function of the Fermi momentum, kf, the 
(K^) condensate, 9, and hb SiS well as the other coefficients in the Lagrangian. First 
we functionally minimize VL with respect to kf{x), and find 



Pb 



r^2 

* I t 



(4.11) 



where /i* is the quasi-particle's chemical potential defined as the value of e* at = kf. 
Eq.( [4.11| ) is equivalent to Eq.( [4.4| ) with k = kf, hence, /is is equivalent to the energy 
of a neutron at the top of the Fermi sea, which is therefore constant over all space. 
Different choices for the chemical potential lead to different sizes of finite solutions, 
i.e. to different numbers of neutrons within it. 

Functionally minimizing Q with respect to 6{x) results in a differential equation 
for the condensate, 9, 

V'(/.^)-^(;^[^n-W,] = 0, (4.12) 
where P„, the pressure of the neutron gas, is given by 



fiBn'n - En 

1„.,- 

3 



CI 



CI 



= (nf ■ kn) + ij:,in^nf - ^(nn? 



2/r ' vr 



In the zero temperature ground state the neutron number density is 

, ^ , X , ^ , fkf d^k k^f 
p„ = {*|„t„|*)=2/ p^ = 5^. 



the scalar density 



PS 



\l>|nn|\l>) = 2^ 

kfPB 



27r2 



d^k m* 



(27r)3 



m* ^ ( He + kf 
—— m ' 



,Pb - kf J 



(4.13) 
(4.14) 



(4.15) 



(4.16) 
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and 



*2 



47r2 



(27r)3 VF + m 





Hb - kfj_ 



■ (4.17) 



Substituting Eqs.([4.5D, ( [4.11| ), ( ^.151 ), and ( [4. 161 ) into the equation for the disper- 
sion relation at the top of the Fermi sea, 



fi*'-k 



f 



m 



*2 



(4.18) 



we get a a transcendental equation of the variables kf and 6 which we solve numerically 
to get kf{9). We can now solve Eq.( [4.12| ) after plugging kf{9) into Eq. (|4.14|) to obtain 

Pn{0) . 

We can view Eq. (|4.12| ) as a one dimensional newtonian equation of motion ^ for 
a particle of unit mass moving in a potential, Ve/r, with the following replacements: 



fJ 



X, 



t, [p^-u^]^v,ff{x). 



(4.19) 



Writing Eq. ( |4.12| ) in radial coordinates, we get (for our spherically symmetrical, time 
independent condensate 9) 



<fUJ) ^ 2d{U 



rdr 



d{Pn - Ue 

dim 



(4.20) 



Notice that the "damping term" decreases as 1/r which means it becomes negli- 
gible for large "times" in the newtonian analogy. The potential, V^ff = Pn — Ue for 
an infinite baryon solution, with Cy = .24, C| = .63, and fis = 900MeV, is shown 
in Fig.l. The potential has two degenerate maxima one at 9q = 1.37 and the other 
at the true vacuum, ^ = 0. In the newtonian analogy a test particle starts from the 
top of the hill at 9o = 1.37 and waits there for a very long time during which the 
damping term becomes effectively zero. The particle then accelerates quickly through 
the valley reaching the top of the other hill, the true vacuum, where it comes to a 
stop. Therefore, the solutions for a large number of baryons have constant 9 and 
density over a large volume and a small spherical surface region where one vacuum 
state evolves rapidly to the other. One can isolate the infinite solutions by imposing 
the following conditions on Vejj, 



VeffiOo) = V,ff{9 = 0) 







d9 



(4.21) 
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where is the value of the condensate at the center of the solution, r = 0. Notice 
that for an infinite solution V^jjIOq) = fisn'^n — En ~ W^,, = and in the constant 
interior of the solution e = e„ + giving the relation, 

liB = -, (4.22) 

Pn 

thus the energy per neutron is /i^. 

The strangeness density of each solution is found as the zero component of the 
strangeness current, S,^i = {Jy)i^i — {Jb)ij., where Jy is the strong hyper charge current 
and Jb is the baryon number current. With our approximations the only non zero 
contribution to 5*0 comes from the meson vector current coupling iTrB'y^[V^, B] term 
in Eq. (|2.4|) , and it is 5*0 = (1 — cos6'o)n"''n so the strangeness per baryon number of 
the strange matter is just (1 — cos 6*0) (see Table 1). 



5. Mapping Out the Solution Space 

The model we are considering so far with only n = 2 terms, which is equivalent to 
the Walecka model for bulk nucleonic matter has only two parameters, Cy and 



C|. Notice that the coefficients C| and Cy are positive definite and the signs of the 
terms are chosen to account for the scalar attraction and vector repulsion observed 
in nuclear interactions. Thus we have four unconstrained parameters, C|, Cy, fiB, 
and 6*0, and two constraining conditions on Vejj, Eq. ([4.21|) . We are left with two 



independent variables that we choose to be /is and Cg. Since the main properties 
of the infinite solutions we are looking for are their baryon density, pn, and their 
binding energy, Ei,, we will show the regions of solutions in the {pn, Ef,) space. This 
is equivalent in the variables we are talking about to a {Oq,pb) space. This is so 
because p„ = k^/Sn'^ depends on 6*0 through the monotonically increasing (see below) 
function, kf{9o), and Ef, = pB — fnN, see Eq. ([4.22|) , inside the infinite solution, i.e. 
in bulk baryonic matter. For a given pB, i.e. a given Ef,, there is a range of allowed 
number densities, that depend on Cy and C|, p„ = Pn{Cy, C|), but Cy also depends 
on C|, Cy = Cy{Cg). Cy cau be shown numerically to increase monotonically with 
C|, while Pn decreases monotonically with C|. So, the allowed range of p„ for each 
Eb corresponds to an allowed range of C| values, < C| < (C|)maa:- It turns out 
{Cg)max can be found in a systematic way, that we pass now to explain. 



Let us return to Eq. ([4.18|) , the transcendental equation whose roots we solve 
numerically for to find kf at a fixed 6. Using Eqs.( [4.1I| ) and ( [4.15| ) we see that 
p* = Pb - (Cv^//^)(fc//37r2), therefore the l.h.s. of Eq.( [4.18| ) does not exphcitly 



depend on 9. Through Eqs.( [4.5| ) and ( [4.16| ) we see that the r.h.s., m* , carries the 



exphcit 9 dependence of Eq.( |4.18| ) through m*. The l.h.s. depends on pB, Cy, and 
kf and the r.h.s. depends on C| and kf. Figs. 2a and 2c are examples of the sides 
of Eq.( |4.18| ) as functions of /c/ for particular C| and Cy values. The l.h.s are shown 
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with black dashed hnes for a given /i^ {fJ^B = 900Me\^), the r.h.s. are shown for 
different 6 values with gray solid lines. The shape of the gray solid curves depends in 
a complicated way on C|. For each 6, the intersections of the corresponding gray line 
with the dashed line gives the solution kf{0). The dashed line moves up and down 
the diagram with without changing shape. For a fixed /is, as ^ increases from zero 
the solution kf{9) starts departing from zero and grows continuously to a maximum 
value where all the gray lines turn over, thus kf{9) is a monotonically increasing 
function. This happens when there is only one intersection (for ^b and 9 fixed). This 
is the case of Fig.2a. However there are cases, such as the one shown in Fig.2c, in 
which there are three intersections (see the three x's). We can see in Fig.2c that as 9 
increases (gray line lowers with respect to the dashed line) the first two intersections 
approach each other, get to coincide at the point where the two intersecting curves 
have the same slope, and disappear, leaving only the third intersection. If we had 
taken the intersection as the physical kf{9), at the 9 value for which the two first 
intersections disappear, after joining in one point, there would be a discreet jump 
in kf from this point to the 3'''^ intersection, at a larger value of fc/, which becomes 
now the only intersection. This jump in kf is unphysical, since the density p„ ~ kj^ 
has to be a continuous function of the condensate, 9, as 9 grows from zero outside 
the bound state to its maximum value in the interior. Thus we impose conditions 
on the l.h.s. and r.h.s. of Eq. ([4.18|) as functions of kf to forbid the values of C| 
(for each yU^) for which triple intersections appear. These are conditions on the slope 
of the intersecting curves (i.e. the l.h.s and the r.h.s.) as functions of kf. Valid 
solutions only occur when at kf = the r.h.s, m*^ (gray line) is lower than the l.h.s., 
{fi*^ — kj) (black dashed line), see Figs. 2a and 2c. This insures that as the condensate 
9 increases (gray line lowers) kf{9) (and the baryon density) increases smoothly from 
zero near the surface of the solution^. The curves cross when they intersect and m*^ 
becomes larger. If the m*^ curve is more concave than the (/i*^ — k'j) one, they cross 
again. This is the case we reject, namely when 

for any value of 9 and kf. Notice that both slopes are negative and we reject the 
case in which the m*^ slope is larger (more negative) than the other. We show the 
difference in slope for a rejected case (the one of Fig.2c) in Fig. 2d. Fig. 2b shows an 
allowed border line case (the one of Fig. 2a) in which the difference in slope would 
become negative (for some values of kf and 9) if C| would be increased. This is the 
important numerical result that allows us to get an upper bound on C|, the fact 
that the difference in slopes (i.e. the l.h.s. of Eq.( |5.1| )) decreases with increasing C|. 
Thus, the border line case for each /i^ fixes {Cg)max- 

^Notice that if instead the black curve is lower than the gray curve at kf — 0, as 9 increases the 
first intersection happens at fc/ ^ and approaches kf ~ 0. This corresponds to the unphysical 
situation of a discontinuous density that starts all of a sudden with a large value near the surface 
and decreases at the center of the solution. 
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Now we can explain the procedure actually followed to map out the solution 
space. First, we choose a /i^ (or equivalently Ef,). Then vary Cy and C| and look 
at Veff (solving Eg. ( [4. 181) to get kf{6)) until the conditions for an infinite solution, 
Eq.( [4.2lD , are fulfilled, namely until the maximum of Veff is at zero. This is not 
difficult to do after noticing that increasing C| (scalar attraction) raises Veff and 
increasing Cy (vector repulsion) lowers Veff. Once a set (C|, Cy) corresponding to an 
infinite solution has been found, we look at the l.h.s. of Eq.( p.lD (i.e. the difference 
of the slopes of both sides of Eg. ( [4.1^ )). If it corresponds to an allowed case (positive 
difference in slope) C| is increased, otherwise C| is decreased, and the shape of Veff has 
to be checked again to obtain another Cy so that the new (C|, Cy) corresponds to an 
infinite solution. One keeps doing this until finding an infinite solution corresponding 
to a border line case (Fig.2b) for the difference in slope. At this point we have 
found {Cg)max for the given /x^, and its corresponding Cy and 9q. This procedure 
determines the {pn)min border of the allowed regions in {Eb, pn) space. The {pn)max 
border corresponds to C| = and the Cy value necessary to satisfy Eq. (|4.21|) , for 
each pb- The allowed regions found are shown in Fig.3 with labels 03 = 1.3 and 
as = 1.5. The parameter 03 appears in m„*, Eq.( |3.3| ). It is measured through the 
nucleon cr-term to be 03 = 1.3 ± 0.2. It is the parameter responsible for the s-wave 
attraction between kaons and nucleons, that yield K-condensation. We see in Fig.3 
that a higher value of 03, i.e. a lower value of the effective nucleon mass, m*, helps 
to find solutions at lower densities. Only values of Cy > are allowed, producing the 
lower boundaries shown in the figure with the Cy = label. Some values of {Cy, C|) 
are also shown. 

Notice that the procedure we follow is only self consistent for pn/ Po < 7, since 
the perturbative expansion in baryon bilinears BTaB is lost at larger densities. The 
rest of the diagram for p^ > 7po in Fig.3 may be indicative but is not believable. In 
any case, we were not expecting solutions below the density necessary for the onset 
of K-condensation, 2-4 po- This is consistent with the densities we find, starting at 
4.5po- At these densities, the fermion bilinear expansion parameter ~ Pn/7po is large 
and higher order terms with more baryons should be considered. 

6. Solutions with Higher Order n = 3 Terms 

Let us describe the effect of taking into account the Ci, C2, and Hi terms in Cchirai, 
Eq.( |3.2|) . These are the only remaining terms of order = 3 (a + c = 3), besides 
the ai,a2,as terms {a + d = 3), already included in our initial lagrangian, because 
they are the lowest order terms providing an s-wave (i.e. derivative or momentum 
independent) K-B couplings, thus without them kaons and baryons would not be 
coupled within bulk baryonic matter. 

The procedure we follow to examine the effect of these higher order terms is to 
add only one term at a time. This increases by one the number of free parameters, 
Pb, Oo, Cy, Cg. However we reduce the problem again to just the old four parameters 
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by fixing the new one. We choose both negative and positive values of 0(1) for the 
single higher order parameter studied, within a range that insures the higher order 
term in the lagrangian is never larger than the Cy and C| terms. Then we proceed 
as before choosing a fiB value and finding (C|)max (and the corresponding {pn)min, 
etc.) as described above. 

We are interested in knowing which higher order terms help to get solutions with 
lower density p„, for a given binding energy. The three baryons Gi and G2 terms 
raise {pn)min instead, as well as the Hi term with positive Hi values, while with 
negative Hi values it can lower {pn)min although very little. An example of the effect 
of these terms is shown for Gi = 0.4 (and = 1.3) in Fig.3. The region of allowed 
solutions in (p„, Ej,) is a wedge because {Cg)max becomes zero at the tip (and the 
{pnjmax bouudary corresponds to C| = 0). Other examples with Gi, G2, and Hi 
positive and negative are given in Table 1. This table shows the values of (C|), (Cy), 
the one higher order parameter chosen to be non-zero, the baryon number density, 
the effective mass m* of the quasi particles and the condensate ^o, for which the 
baryon number density is minimum (corresponding to the maximum values of (C|) 
and (Cy) for which solutions exist, after fixing the higher order parameter to the 
largest positive and negative physically acceptable values) for a fixed binding energy 
El, = —39MeV and parameter 03 = 1.3. While Gi does not introduce any major 
change in the self consistent procedure described above used to find solutions, the 



G2 and Hi terms turn the p* algebraic Eq.( 4.11 ) into a transcendental equation, (see 
Eq.(A.l) in the appendix with either ifi = Jy = or G2 = /y = 0). The appendix 
gives the complete equations with all the higher order terms included in this paper for 
the effective neutron chemical potential, effective neutron mass, pressure and energy 
density that we label with the sub index ''total" ■ Notice that the Hi term in Eq.(|3.2|) 
contains a derivative of the nucleon field so it modifies the momentum of the neutrons. 
Adding only Hi to the Cy and C| parameters, Dirac's equation becomes 



H -> H -> 

7°e* - (1 + -j^inn))^ ■ k - {m* - -j^{n^ ■ kn)) 



n = 0. (6.1) 

where e* and m* are those found in Eqs. ([4.4| ) and ([4.5|) . One can arrive at a standard 
free particle dispersion relation by dividing Eq.( |6.1| ) by the factor multiplying 7 -A; and 
redefining the effective energy and mass. This is the origin of the (1 + {Hi/ f^A)nn) 
denominator in /x^oj,^^ (i.e. e*TOTAL at the top of the Fermi sea, where k = kj) and 
'"^TOTAL given in the appendix Eqs.(A.l) and (A. 2). 

The importance of having to solve a transcendental equation for p*{kf) when G2 
or Hi are non-zero is that many of their solutions are incompatible with the existence 
of infinite solutions. The reason is that solutions for p*{kf) can not be found for large 
ranges of G2 or Hi for p*'^{kf) < kj- — A, where A depends on G2 or Hi and Cy. 
Since /i*^ — kj = m*^ (where e* and p* are e^oTAL ^71*^^4^ given in the appendix 
with only C|, Cy, and G2 or Hi respectively non-zero) this means that the effective 
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mass has a non-zero minimum m*^ > A = m^^^ > 0. This forbids solutions that 
would require m*^ < A in the bulk matter. Since A increases with increasing values 
(positive and negative) of G2 and Hi, only small values of these parameters lead to 
solutions, and even then {pn)min is not improved. 

From this analysis we conclude that pure nuclear interaction terms do not help to 
lower {pn)min for a given binding energy. Then, one is lead to try higher order terms 
dependent on the condensate. As we will see they may actually help. 



7. Solutions with Higher Order Condensate Dependent Terms 

The first terms of this type are of order n = a + d = A. They are obtained by 
squaring the trace in the 2'^'^ term of Ceff in Eq. (|2.4|) , by multiplying the ai, 02, and 
03 terms in Eq.(pl^ by BB/f^A, by multiplying the Cy and C| terms in Eq.(|2^) 
by Tr[A1E + h.c] and by writing similar terms with alternative ways of taking the 
traces. They all give only the following three terms within bulk matter with only 
neutrons in a condensate, 

~ J> \ ~ cos 6)nnnn ms{l — cos 6)n''nn''n 

-Jf^m^il-cos'O). (7.1) 

The terms Is and Jy effectively add a 9 dependence to the Cs"^ and Cy^ terms, 
respectively. While ly, that modifies p* (see Eq.(A.l)), does not help to lower 
positive values of Is, that modify m* (see Eq.(A.2)), do help (see Table 1 and notice 
that {pn)min = 5.3po for Is = 3 instead of 6.5po for Is = 0). The J term helps in 
lowering p„, and it is experimentally determined, because it corrects the kaon mass 
and other terms in the potential energy, Ue (see Eq.(A.5)). We must be careful since 
now mj^ = rrisA + ... because we include higher order terms in m^. Let us now call 
friK the value of the K-mass at first order in m^, namely mj^ = rUsA. Therefore we 
change the notation which we used in Eq. ( ^.2[ ) for the term /^m|'(l — cos 6') into 
flm\{l — cos 9), and we write U9 as 

Ue = f^ml{l-cos9) + J^m%{l-cos''9) + ---. (7.2) 



In order to evaluate m^, the 1** order value of the kaon mass, we use the Gell- 

4 = 



Mann-Okubo relation, that holds only at first order in the quark masses, m'j 



l/4(3m^ + m^), where we take and to be the physical masses of the vr and 
7] mesons respectively. We find = m^A ~ (ASOMeV)"^. We take this as a rough 
estimate of the minimum value of mx-f\ 

^Taking ttik instead of the physical tuk mass to be equal to m^A amounts to a redefinition of 
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Expanding cos^ 6* and leaving only the terms quadratic in 9, = 2{K^)'^ / f^, 
Eq.([7.2|) becomes 

\\mKJ A^mj^ J 
= m^(ir°)'(0.94 + .43J + ■■■) + ■■■, (7.3) 

where rriK is the physical kaon mass, ttik = A98MeV. Since the parenthesis must be 
1, we get an upper bound on J, J < 0.14. 

In order to examine the most favorable case, the lowest p„ for a given Et, we take, 
besides the largest reasonable value for 03 i.e. 03 = 1.5, J = 0.14 and Is as large 
as possible, while keeping the Is term smaller than the C| term in agreement with 
a perturbative expansion. We show the results in Fig. 3. Because the largest value 
of ^0 encountered in the corresponding solutions is 9q ~ 1.5, the value Is = 4.6C| is 
the largest insuring the Is term is never larger than the C| term. These solutions, 
those with the lowest p„ we find for every binding energy, are our main result. We 
believe the perturbative chiral lagrangian cannot reasonably do any better. Notice 
that solutions with p„ < 4po have at most only 50MeV of binding energy, and those 
with pn < 5po do not have more than Eh — 70MeV. If we were to accept densities 
closer to 7po we could hardly get to Ej, ~ 150MeV^. These binding energies may be 
too small to allow for the formation and persistence of this type of strange baryonic 
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8. Comparison with the LNT Model 

Since much lower densities, even 1.4po, and larger binding energies, up to Eb — 
300MeV, have been found by LNT [|l^, we want now to discuss the origin of this 
difference. 

LNT take the lagrangian in Eq. (|2.4|) but without the four and six fermion terms, 
add to it the lagrangian of the Walecka model [jlT|, a renormalizable lagrangian intro- 
ducing two fictitious fields (p and coupled to baryons, and add the following term 
coupling (f) to the pseudo-Goldstone bosons 

^f^A[TrM{J: - 1) + h.c.]{h(j) + h2(f' + ■■■), (8.1) 

where hi and 62 are arbitrary constants and (f) (and V^) are taken to be chiral singlets. 
At a first glance this coupling seems to violate the principle of using QCD compatible 
M. dependent terms, obtained by promoting to a field that transforms under a 
chiral transformation. This trick apparently would forbid terms proportional to TrM. 



the strange quark mass TOj. This, however, would not effect our previous results because they do 
not depend explicitly on m^, but rather on such combinations as m^as, mgls, etc, and the change 
in rris can be compensated by small changes in the accompanying parameters. 
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in the lagrangian. It is true that chiral lagrangians do not contain (j) fields and the 
rules to construct them do only apply to mesons and baryons. However, notice that 
when the (f) field is heavy, through (f) exchange at tree level the hi coupling in Eq.( |8.1| ) 
generates an effective chiral lagrangian coupling proportional to [TrA1(S — l) + ./i.c.]^, 
which contains a TrAi term. However, at a second glance one can see that one can 
start from the coupling dictated by the trick, namely ~ Tr[A1E + /i.c.](6i0+620^ + ' ■ ■) 
and by shifting the field to the minimum 0o of the potential for E = 1, namely 
replacing by (f)+4>o, one obtains the lagrangian used by LNT, after a few redefinitions 
of constants. 

LNT use the Walecka model to account for the existence of nuclei, not included in 
the chiral lagrangian they study, which does not have four fermion terms. However 
it has been shown that precisely these terms are equivalent to the Walecka model in 
bulk nucleonic matter |]l3l . Thus it is not necessary to go out of chiral lagrangians to 
account for normal nuclear matter. Besides, the additional couplings in Eq.( p.4D of 
the Walecka fictitious fields with mesons are arbitrary. 

We here choose instead to use solely perturbative chiral lagrangians. We believe 
this is the main difference. Although the LNT lagrangian is chirally symmetric it 
does not seem to be obtainable from a perturbative chiral lagrangian. We will now 
compare their solution with the chiral expansion we used so far. In order to do it, we 
need to eliminate from the LNT lagrangian for bulk matter, by using the equation 
of motion for constant, 

, ^ g^nn - f^mj^bijl - cos 9) 
The LNT lagrangian becomes 

c = ^f:d,ed>^d-f^ml{i-cose) 



+n{i<f) — mN)n + atms{l — cos9)nn 



/ t„^2 I l [gcf>nn- f^mlbijl- cose)] 
2ml^ ' ^2 ml + 2f^mlh{l-cose) ' ^ ' 

where at = + 02 + 01 = 2.32 ± 0.4, 171^ = 939MeV is the mass of the free neutron 
(see Eq. (|3.3|) ) and g^i,, m^j,, gy, my are respectively the couplings to nucleons and the 
masses of the scalar and vector Walecka fields. 

Using the procedure described above (the same used by LNT) one finds the LNT 
effective chemical potential and effective mass of the quasi-particles in bulk matter, 

f^* = f^B-^n^n (8.4) 
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I, .X glnn + g^flm\hi{l - cosO) 
m = TUN- atms{l - cose) 2,0/2 2 u n o^' ^^■^f 

These formulas reproduce the results of LNT. 

The main feature of the LNT model that allows to obtain very low p„ and large 
Eb is that the value of at is at = 8.2 ± 0.2 instead of at = 2.3 ± 0.2. We have already 
seen that large values of 03 help greatly in reducing the nuclear density of solutions. 
This large value of 03 is obtained by LNT, because the exchange of the particle 
contributes to pion-kaon scattering so that now 

as- ^'^^"' =L3±0.2, (8.6) 
2m^ 

which for bi = 20GeV^~^, the value chosen by LNT, gives 03 = 4.2 ± 0.2. Note that if 
the last terms of Eqs.(|8.3|) and ( p. 51 ) are expanded in powers of (1 — cos6')'' and the 
term in (1 — cos^) is summed to the at term we obtain precisely Eq.( ^.6|) (one needs 
to refer to the original lagrangian to see that only a^ and not ai and 02 are corrected). 
Let us write the terms in Eq. (|8.3| ) that depend solely on 6 (not on n). These will be 
enough for our argument. 

Ue - f.mAl - cose) - ^ 2f^mlb,{l - cosey ^^''^ 

To rewrite these terms in a form similar to chiral perturbation theory we use the 
Taylor expansion 



;i + f3^i _ cose))-' = (1 + P)-' I 1 + E (^Y 



COS-' e 



where /3 = 2f^b2{mK/Tn^)'^ ~ 1.66, to expand the denominator of Eq.( p.7| ) which 
becomes, 

= -f.<(i-c^P^]cose 



where 



1 f>W 

2{1 + P)ml 



a = - //^ g/^ -0.62, (8.10) 



once the constant term has been dropped. 
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In chiral perturbation theory the terms of IAq of power j in cosO are proportional 
to {rris/Ay = 0.24-5', go naively we would expect each term to be 0.24 of the previous 
term. Moreover Eg. ( |7.3| ) shows that the first order contribution to the kaon mass 
accounts for 0.94 of the total mass. Let us isolate the kaon mass term in Eq. (|8.9|) by 
expanding cos^ to order 6*^ = 2{K^)'^ / , 



2 + p a 
a- 




+ ■■■ (8.11) 



The terms in the parenthesis of Eg. ( |8.11| ) sum to unity as we expect from Eg.( |8.7| ) 
(where the only 6"^ term is m\{K^)'^). However, writing out the first few terms 
exphcitly in Eq. (|8ll|) , 

Uq^'^ ~ m^(ir°)^(0.15 + 0.18 + 0.16 + 0.14 + ■■■) + ■•• (8.12) 

we see that the series converges very slowly in comparison to the perturbative expan- 
sion in chiral lagrangians shown in Eg. ( [7.3|) . 

Another way of showing the same difference between the LNT models and ours 
is by obtaining the LNT coefficient hi in terms of the chiral expansion coefficient J. 
This is most easily done by expanding Ue in powers of (1 — cos 6*) instead of powers 
of cos 6 

Ue = {flml + 2J/^^^)(1 - cos9) - J fl^{l - cosOf + ■ ■ • . (8.13) 

Writing lAe in this way it becomes obvious that the higher order term for J positive 
lowers the potential energy of the condensate. Notice the factor multiplying (1 — cos 9) 
in the first term, is f^mj^ as shown in Eg.( |7.3| ). 

Referring back to Eg.( |8.7| ), we also expand C in terms of (1 — cos 6') (only for 
6 < 1.16, so that 2/2771^62(1 — cos6)/m?^ < 1) and obtain 

U,^^^ = />2^(l-cos^)-4?^(l-cos^)^ 



2mi 



m2 ^2 



Comparing the coefficients of the terms proportional to (1 — cos^)^ in Egs.( 8.lS ) 
and ( ^.14| ), from the bound J < 0.14 we obtained above we get bi < 2.7GeV^ , 
which is much smaller than 20GeV~^, the value used by LNT. Thus, even if the LNT 
lagrangian is chirally symmetric it does not seem to be obtainable from a perturbative 
chiral lagrangian. 
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9. Conclusions 



The main results of this paper are given in Fig. 3 where we show the lowest baryonic 
number densities pn corresponding to infinite solutions with a given binding energy Ei,, 
that we obtain with a chiral effective lagrangian, as explained in section 7. Actually 
only the region of p„ < 7po is consistent with our approach. In section 5 we found that 
the lowest order SU(3)lXSU(3)j^ chiral lagrangian, including four fermion terms, could 
not yield solutions with densities smaller than 4.5po (see Fig. 3) and at these densities 
the fermion bilinear expansion parameter BTaB/ f^A ~ Pn/7po is large and higher 
order terms with more baryons should be considered. We analyze them in section 6 
and find that they do not help to lower {pn)min for a given binding energy. In section 
7 we show how much higher order condensate dependent terms help, when their 
coefficients are chosen to respect the perturbative expansion. We find that solutions 
with pn < 4po(7po) have at most only 50MeV{150MeV) of binding energy. We believe 
the perturbative chiral lagrangian cannot reasonably do any better. The densities we 
obtain are entirely compatible with the onset of a K-condensate at densities 2-4 pQ. 

The binding energies we find may be too low to allow for the formation and 
persistence of this type of strange baryonic matter in the Universe at temperatures 
of order lOOMeV [|T5|. However, we have used only a f^^-condensate and neutrons, 
the lightest baryons in such a condensate. We may not exclude that a less restrictive 
ansatz may provide solutions with larger binding energies (a factor of 3 or 4 may 
suffice lO) that may appear and persist in the early Universe. 
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Appendix A 



These are complete equations with all higher order terms included in this paper for 
the effective neutron chemical potential, effective neutron mass, pressure and energy 
density, that we label with the sub index "total"- See Eqs.( [4.15D , (^4.16|) and ( |4.17|) 
for (n^n), {nn) and (n7 ■ /cn), respectively. 



/i. 



TOTAL 



■ PA 



'1 — cos 9){n^n) 



1 + ^(nn, 



m 



N 



Q2 
fn 



fA 

' PA 



nn) 



(n^n) 



fA 

cos 6){nn 



l + ^inn) 



n -1 



(A.l) 



(A.2) 



1 ,^ AHi _ 



Cl 



- (1 + -TTT-nn) (n7 ■ kn) + — — - — — 



PA 



2/2 2/2 3/,4A 



[nn] 



2G2 

'PA 



{nn){n'^ nY H 7^(1 ^ cos^)(n^r2) 

J TV 



VA 



'1 — cos 6){nn) 



(A.3) 



2^1 ,3 



+77T7v(^^)' + ^^^(1 - cose)(ntn) + ^^(1 - COS e){nn)+Ue 



3fA 



VA 



VA 



(A.4) 



Ue = msAf:{l-cose) + Jpmiil-cos 



(A.5) 
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{Cg^jfiax 




n > 3 


{Pn)min/ PO 


m*{MeV) 


Oo 


0.62 


0.24 






6.6 


220 


1.4 


0.00 


0.17 




0.65 


8.1 


180 


1.7 


0.77 


0.21 


-1.00 


7.1 


200 


1.6 


0.61 


0.27 


G2 


0.04 


6.6 


230 


1.3 


0.52 


0.01 


-0.25 


6.9 


210 


2.0 


0.58 


0.05 


Hi 


0.50 


7.8 


190 


1.9 


0.64 


0.28 


-0.10 


6.3 


240 


1.3 


0.34 


0.02 


Iv 


0.80 


7.3 


342 


1.3 


0.89 


0.48 


-0.90 


7.6 


73 


1.8 


0.72 


0.36 


Is 


3.00 


5.3 


210 


1.2 


0.58 


0.19 


-1.00 


7.4 


230 


1.6 


0.66 


0.26 


J 


0.14 


6.4 


210 


1.4 



Table 1: This table shows the values of (C|), (Cy), the one higher order parameter 
chosen to be non-zero, the baryon number density, the effective mass m* of the quasi 
particles and the condensate Oq, with which the baryon number density is minimum 
for a fixed binding energy Eh = — 39Mel^ and parameter = 1.3. This corresponds 
to the maximum values of (C|) and (Cy) for which solutions exist, after fixing the 
higher order parameter to the largest positive and negative physically acceptable 
values. The first entry is the solution in Figs.l and 2. The strangeness/baryon for 
these solutions is (1 — cos^o) which lies in the range 0.64 - 1.4. 
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Figure 1: An infinite solution with binding energy = —39MeV, number density 
Pn = 6.6poj and parameters Cy = .24 and C| = .62 (and no higher order terms). This 
is a minimum density solution where p = Pmin, Cy = {Cy)maxi and C| = {C'g)max- It 
corresponds to the first entry in the Table 1. 
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Figure 2: (a) The l.h.s. and r.h.s of Eg. ( [4.18D for the solution shown in Fig.l are 
plotted here. The dark dashed line corresponds to the l.h.s of Eq. (|4.18| ), and the 
gray lines to the r.h.s side for the indicated values of 9. The intersection of these 
curves determines the numerical solution for kf{6). (b) This is a plot of the function 
2m*{dm*/dkf)-2^*{dij,*/dkf)-2kf, which is the l.h.s. of Eq.(^, for 6 = 0.4. This 
shows why the solution in Fig.l has minimum density, since for any higher values 
of Cy and C| the function shown here would have negative values and, hence, the 
solution would be unphysical. (c),(d) Show the same as (a) and (b) respectively, 
but for an unphysical solution with the same binding energy, with parameters Cy = 
.27 and C| 



.77, larger than {Cy)max = -24 and (C|)max = -62 corresponding to the 



case in (a) and (b). 
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P"/Po 

Figure 3: Here we show the regions of physical solutions in {Ef,, pn) space correspond- 
ing to the indicated parameters. We indicate only the values of as and the higher 
order parameters taken to be non-zero. We then find the values of C| and Cy at the 
left and right boundaries of each band with the procedure described in the text. C| 
and Cy are maximum at the left boundaries and C| = at the right ones. The lower 
boundary corresponds to Cy = 0. The values of (C|,Cy) are shown at some points. 
We have chosen the most favorable values of and the higher order parameters to 
obtain the lowest baryonic number densities p„ for infinite solutions with a given 
binding energy E^, using a chiral effective lagrangian. Actually only the region of 
Pn < 7po is consistent with our approach (the perturbative expansion is lost at higher 
densities), so the rest of the figure may be indicative but is not believable. 



26 



